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Abstract. Results on direct photon measurements from the PHENIX experi- 
ment at RHIC are presented. The direct photon yields for >6GeV/c as a 
function of centrality in Au-Au collisions at % /J/v r /v=200 GeV are found to be 
consistent with NLO pQCD calculation scaled by the number of binary col- 
lisions. The results suggest that the photons observed are emitted from the 
| initial stage of hard scattering. Comparisons with several theoretical calcula- 

tions are also presented. 
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O \ 1. Introduction 



The Relativistic Heavy Ion Collider (RHIC) at BNL is expected to produce a hot 
and dense matter, whose energy density is sufficiently high to induce a phase tran- 
sition to the so-called Quark Gluon Plasma (QGP) [[I]. RHIC has been operated 
for five years since the first collision of Au+Au in 2000, and produced many intrigu- 
ing results that have never appeared in lower center-of-mass energy systems [ ■ 
Especially the energy loss of hard-scattered parton in the hot and dense medium 
that results in "jet quenching" has been of great interest since its first discovery by 
the PHENIX experiment [|3J. The hard scattered probe now become an useful tool 
to explore system dynamics because its cross-section is large at RHIC energy, and 
can be accurately calculated by perturbative QCD (pQCD). The results obtained 
so far including jet quenching suggest that a hot dense medium has been created. 

Photons are excellent electromagnetic probes for extracting the thermodynam- 
ical information of where they are emitted because they do not interact strongly 
with medium once produced. Thus, they are expected to provide hints to answer 
questions of whether the jet quenching is due to initial or final state effect, and if the 
medium turned into the QGP. The photon production processes are shown in Fig.^ 
The Compton scattering of quarks and gluons and the annihilation of quarks and 
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Fig. 1. Photon production in (a) leading order process, and (b) next-to-leading 
order process. 

anti-quarks are leading order processes, and the next leading process is dominated 
by bremsstrahlung. There is also a prediction of a jet-photon conversion process, 
which occurs if QGP is formed, by a secondary interaction of a hard scattered par- 
ton with thermal partons in the medium Fig. shows a theoretical prediction 
of photons over apr range The calculation predicts that a photon contribution 



> 

CD 

a 



Central Au+Au (s =200AGeV) 

<N dl >-800 
|y|<0 35 



Central Pb-Pb s ~=17.3A(3eV 



Hadron Gas 
QGP (T =205MeV) 
pQCD (Ak t 2 =0.2) 
Sum 




WA98 Data 

Hadron Gas \ P 10 til T-L V^St, 

QGP (T=370MeV) \ S ** Vv 5 ***. 2.35<y<2.95 

initial pQCD (pp) 
sum 

1 2 

q, [GeV] 

Fig. 2. Theoretical calculation on various 1 ,„ 3 4 

° q, [GeV] 

photon contributions for central Au+Au „. _ ,, 7 . AO , , , . ,, 

„. . , Fl g- 3 - WA98 results and their thc- 

colhsions at ,/i]v]v=200GeV. ,. , . , 

v oretical interpretation. 

from the QGP state is predominant in the pr range of 1< pr <3 GeV/c. The signal 
can be seen after subtracting photons decaying from known hadronic sources. The 
typical signal to background ratio is ~10% [E]. For px >3GeV/c, the signal is 
dominated by a contribution from initial hard scattering, andp-r <1 GeV, the signal 
is from hadron gas through processes of irir(p) — ► 7/5(71"), TtK* — > and etc.. 

Direct photon measurements in heavy ion collisions were first carried out by the 
WA80 experiment in 200AGeV S+Au collisions, where only the upper limits are 
obtained [[7]. The WA98 experiments observed a significant excess in the pr region 
of 1< p T <3GeV/c in 158AGeV Pb+Pb collisions [|U E], and a recent theoretical 
calculation interpreted it as a consequence of either kr broadening or high initial 
temperature [E3j as shown in Fig.|3 Since there is no direct photon result for p+Pb 
collisions, one can not evaluate the effect of kr broadening in Pb+Pb collisions. 
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The E706 Collaboration at Fcrmilab measured direct photons in 7r+Be, p+Be at 
515 and 800GeV/c, and showed that the data deviate from NLO pQCD calculation 
without including additional kx smearing according to the incident particle and 
target species [lll|. The result suggests that the Cronin-like enhancement can be 
explained by the initial fcr broadening. The STAR experiment at RHIC measured 
inclusive photons in the pt region of 0.5< pt <2.5GeV/c, and compared it to the 
possible decay 7 contributions from ir° [ ■ Since the systematic error is order 
of 40-50% in the measurement, it is not possible to conclude whether there is a 
significant excess or not. 

In this paper, the latest results on direct photon measurement from the PHENIX 
experiment are shown, and the source of the photon contribution is discussed. 

2. PHENIX detector 

Figure |3| shows the PHENIX detector at RHIC. Both central arms (East and West) 




West 



Fig. 4. The PHENIX detector. 

of the detector include a Multiplicity Vertex Detector (MVD), Drift Chamber (DC), 
Pad Chamber (PC), Ring Imaging Cerenkov Counter (RICH) and Electro-Magnetic 
Calorimeter (EMCal), and cover a rapidity range of |ry| <0.35 and a quarter azimuth. 
The East arm also has a Time Of Flight (TOF) detector and Time Expansion Cham- 
ber (TEC). The detailed description of the detector is given in the literature [ I13j . 
In this analysis, EMCal was used for measuring the energy of photons. The EMCal 
consists of six lead-scintillator sandwich type calorimeters (PbSc) and two lead- 
glass homogeneous type calorimeters (PbGl). The DC and PC were used to track 
charged particles to estimate charged hadron contamination in EMCal clusters. 

The analysis of direct photons require determining background photons decay- 
ing from known hadronic sources such as ir° and r\. PHENIX has measured the 
transverse momentum spectra of 7r° up to 13 GeV in Au+Au collisions [ ^] ITS) 
and p+p collision [^] at v / s]v]v=200 GeV. The momentum spectra of r\ and other 
hadronic sources were estimated by using a fit to ir° pt spectrum and applying mj 
scaling: pr — > {pt 2 — M^ 2 + M% ) l l 2 . r\ has been measured but only tj/tt , which 
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Fig. 5. 77 to 7T° ratio for p+p, d+Au and Au+Au collisions at GeV. 
The average value of 0.45+0.05 is obtained. 

is 0.45+0.05 as seen in Fig. [3] was used to determine the normalization factor of 
the TOT-scaled spectra relative to 7r°. Using above fits, the number and spectra of 
photons decaying from ir° and r\ are obtained. 

The inclusive photon spectra were reconstructed by applying several photon ID 
cuts on the measured cluster energy distributions, and correcting for hadron contam- 
ination and PID efficiency. In order to identify photons, a PID likelihood function 
is calculated from several quantities measured by EMCal such as shower shape or 
ratio of energies among towers, and apply a threshold. It resulted in a significant 
reduction of hadronic clusters in the sample. Since the ratio of 7 to 7r° cancels their 
common systematic errors, the excess of the measured photon over the estimated 
background photon is evaluated in terms of i? 7 = (7 /Tr°) measured /(~f / it ) estimated 
(double ratio). The final systematic error on the double ratio is ~ 12-16%, and 
~ 15-20% on direct photon spectra in Au+Au collisions [I17|. 



3. Results 

Fig. shows the direct photon spectra in p+p collisions together with the NLO 
pQCD calculation with three different cut-off scales The statistical error is 

shown as error bars, and the systematic error as boxes. The data are well described 
by the calculation. Fig. shows a compilation of comparisons of data to NLO 
pQCD calculations in p+p and p+A collisions [EH- As are most of the data points 
in the figure, the PHENIX p+p data points shown in red are consistent with pQCD 
expectation within a factor of 2. The result shows that the NLO pQCD calculation 
describes well the yield of direct photons from the initial hard scattering process 
in p+p collisions at -^5=200 GeV. Fig. [8] shows the 7/71" double ratio over various 
centralities in Au+Au collisions together with NLO pQCD calculation scaled by the 
number of binary collisions in red lines. The shaded bands on the data points show 
the systematic errors, and the bars are the statistical errors. The shaded bands 
around the lines show the uncertainty of the pQCD calculation. The magnitude 
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Fig. 7. Compilation of comparisons of data to 
NLO pQCD calculations. PHENIX data are 
shown as red points. 



of the excess at and above px ^4GeV/c increases with increasing centrality of 
collisions, and is consistent with the calculation. Fig. [5] shows the direct photon 



spectra extracted as jdv 
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Fig. 8. 7/-7T double ratio in Au+Au colli- 
sions. 
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Fig. 9. Direct photon spectra in 
Au+Au collisions over centralitics. 



pQCD calculation scaled by the number of binary collisions. It is clearly seen that 
all nine centralities including minimum bias are well described by the calculation. 
Considering the fact that the calculation in Fig. [3] took the suppression of ir° into 
account which results in less background photons in higher centralities, and that 
the Fig. |5] shows the good agreement between the data and the calculation, we 
conclude that the suppression of 7r° is not attributed to the initial hard scattering 
process, while the direct photon is. Fig. ITUl shows nuclear modification factor Raa 
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as a function of centrality, represented by N part , which is defined by: 

R (l/N AA evt )d 2 N AA /dp T dy 
AA < N coU > /a™ 1 ■ d s a pp /d PT dy ' 

The dotted lines show the uncertainty on the number of binary collisions, and the 
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Fig. 10. Nuclear Modification Factor (R AA ) for integrated yield of pt >6GeV for 
7T° (open circles) and direct photons (solid circles). 

boxes shows the uncertainty of the NLO pQCD calculation. The sold circle points 
are for ir° [ 114) , which clearly shows a suppression with increasing centrality. In 
case of the direct photons shown in open circles, the yield follows the number of 
binary collisions scaling. It once again proved that the suppression of ir° is final 
state effect. 



4. Discussion 

There have been several theoretical calculations on photons from various contribu- 
tions. Fig.llllshows the comparison of 0-10 % central PHENIX data with three NLO 
pQCD calculations [ 1201 121) . two of which have same formulation with different kx 
smearing values, and with the calculation including jet-photon conversion process [ 
Ej. The NLO pQCD calculations with different k T smearing values show a difference 
in the yields below px ~ 6GeV/c, where the data points have large error and can 
not distinguish between models. The difference of yield decreases as pr increases 
because of less steep slope at higher pj>. The model including jet-photon conversion 
describes the data very well for pt >4GeV/c. Fig. 1121 shows the comparison of the 
same data with the model with jet-photon conversion to its components. It can be 
seen that the jet-photon conversion process dominates a px region of 4-7GeV/c, 
therefore, it should dominate the integrated yields for pt >6GeV/c. Accordingly, 
this yield should follow the number binary collisions scaling. Thus, the contribu- 
tion from jet-photon conversion process should scale with the number of collisions 
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PHENIX data (0-10% Au+Au 200 GeV) 
All Photons (Fries el al.) 
Jet-Photons (Fries et al.) 
LO pQCD Direct-Photons (Fries et al.) 
Bremsstrahlung-Photons (Fries et al.) 
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Fig. 11. Comparison of 0-10% central Fig. 12. Comparison of the same data 
PHENIX data with NLO pQCD calcula- with the model with jet-photon conver- 
tions and jet-photon conversion included S10n to its components, 
model. 

as well. However, the process assumes the existence of a hot and dense medium 
or QGP, which is not the case in peripheral collisions. The model does not give 
consistent description of collisions for all centralities. The data are also compared 
with calculation on Fig. and shown on Fig. Within the current error, the 




p T (GeV/c) 

Fig. 13. Comparison of the same data with the calculation shown in Fig. 

calculation is consistent with the data up to 4GeV/c, and starts deviating from 
the data above 4GeV/c. However, It is possible that the fact that pQCD calcula- 
tion in the plot does not include an additional fcy smearing can in part explain the 
deviation. 
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5. Conclusion 



Results on direct photon measurements from the PHENIX experiment at RHIC 
are presented. The direct photon yields for £>T>6GeV/c in Au-Au collisions at 
Y / s]vjv=200 GeV are found to be consistent with NLO pQCD calculation scaled by 
the number of binary collisions. The results suggest that the photons observed 
are emitted from the initial stage of hard scattering. Comparisons with several 
theoretical calculations are also presented. 

The direct photon analysis in d+Au collisions and high statistics Au+Au col- 
lisions with the goal to reduce systematics error and disentangle the contributions 
from QGP state is now ongoing. 
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